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Q ' The two-dimensional flow of viscous incompressible fluid in the domain between two concentric circles 
O , is investigated numerically. To solve the problem, the low-order Galerkin models are used. When the inner 
circle rotates fast enough, two axially asymmetric flow regimes are observed. Both regimes are the stationary 
flows precessing in azimuthal direction. First flow represents the region of concentrated vorticity. Another 
flow is the jet-like structure similar to one discovered earlier in experiments [HIS]. 
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Introduction 



i The interest to the problem considered is stimulated by the strange flow regimes induced in 



the thin fluid layer by the rotating cylinder, experimentally discovered by V. A. Vladimirov [H E] 
in 1994. The flow in a thin annular gap restricted in axial direction by rigid walls was considered. 

^ . When the inner cylinder rotates fast enough, axially asymmetric, slow-precessing, stable jet-like 

^ '. pattern appears. These jets represent the sectors with the strong radial outflow which precesses in 

O ■ azimuthal direction. 

^1 The flow regimes observed in experiments are not deeply investigated and the mechanism of 

r-| ■ their onset is not clear. The flow observed is not the Moffatt vortex that appears in the domain 

Ph. with corners [3]. It is not the Hammel flow with the azimuthal jets, investigated and classified 

'. by M. A. Goldshtik, V. N. Shtern et al. (see for example [H [5l [6] where the fiow regimes without 

^ ; the inner cylinder but with the source are considered). It is not the spiral structure described in 

\Q ■ works by M. V. Nezlin and E. N. Snezhkin [71 [8] that appears in the thin layer of rotating fiuid 

^ ■ with a free surface with differentially rotating parabolic-shape bottom. We notice that in spite of 

^ . visual similarity with regimes from [7l Ej, flows in [U [2] are totally different. In [71 [8] the reason 

' for the spiral structure onset is the parabolic bottom proflle that models the gravity fleld in the 

O : radial direction, so that the phenomenon could be described with the use of shallow water theory. 

o\ 
o 



In other words the jet-like structures in [H |8j are similar to waves on the free surface. Here it is 
worth mentioning works by V. Yu. Liapidevsky (see for example [9]) on the two-dimensional vortical 
shallow water flows in a gap between two rigid boundaries. 

There is a number of works on the similar problems. These works are mainly devoted to study the 
5^ ', axially symmetric regimes in the Couette-Taylor flow between short cylinders. These problems were 
intensively investigated numerically after T. B. Benjamin and T. Mullin pLOj. The axially symmetric 
flows in the annular gap with comparable radial and axial sizes, governed by the two-dimensional 
Navier-Stokes equations were studied with the help of the finite-difference methods (particularly 
the methods of markers and cells, and their modifications) and of the Galerkin method (see for 
example [iQl [HI [121 [131 [H]). The base axially symmetric fiow regime in the domain with two rigid 
end-walls represents two Taylor vortices. The bifurcation of this fiow to the regime with one major 
vortex and one minor vortex is described. 

In works [15], [TH [TTl [H] the asymptotic model describing the base axially symmetric regime 
in the thin in axial direction annular domain between two cylinders is constructed and compared 
with experiment. This asymptotic model describes two Taylor vortices in the case of rigid end-walls 
and one Taylor vortex in the case when one of the end-walls is the non-deformable free surface. 
Particularly, the analytic formula has been obtained for the azimuthal velocity (it was earlier defined 
numerically in [HI [121 [13 [H])- The azimuthal velocity calculated with the help of this analytic 
formula is in a good agreement with the experimental data [HI [HJ [171 [18] . 
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In the work presented we attempt to model some of the experimental results from [U !2] with 
the use of the low-order Galerkin model, applied to the viscous flow in the two-dimensional annular 
strip. When the base flow similar to the Couette flow loses stability for the low-order model, 
the stationary asymmetric flow (the soliton-like vortex) appears. This vortex slowly precesses in 
the azimuthal direction. Furthermore, numerical experiments with the low-order models allowed 
to discover another axially asymmetric stationary flow regime: slow precessing in the azimuthal 
direction jet-like structure. This axially asymmetric solution (apart from the first one) does not 
bifurcate from the base flow but appears 'from the sky-blue'. The jet-like regimes appear to 
be similar to ones observed in experiments presented in [HE]- It is rather surprising because, 
according to [H \16\, \T7\ [18], the jet-like flows observed in the narrow (in axial direction) 

gap are substantially three-dimensional. Another surprise is that the Reynolds numbers for which 
asymmetric structures appear are of the same order as in the experiments. Comparison of the 
other parameters is impossible since the numerically solved problem is two-dimensional while the 
experimental facility is essentially three-dimensional. 

The certain uneasiness is caused by dependence of the results of numerical modeling on the 
number of radial basis functions for Galerkin model. With the increase of the number of this 
functions, increase of the critical Reynolds numbers is observed (note that the results weakly depend 
on the number of azimuthal basis functions). In other words, results obtained (existence of the 
precessing stationary flows in the form of the soliton-like jet or the soliton-like vortex) could be 
specific for the low-order Galerkin models only, while the base flow could be stable for the high- 
order models. However, good qualitative agreement with the experiments described in [H [2] make 
us hope that the low-order models are useful to improve the understanding of the mechanism of 
jet-like structures onset in the gaps narrow in the axial direction. 



1 Basic equations 

Let the annular strip domain bounded by the circles of radii ri n r2 (ri < r2) is fllled by the 
viscous incompressible fluid. The inner boundary r = ri is rotating with the constant azimuthal 
velocity fl, and the outer boundary r = r2 is fixed. The Navier-Stokes equations in the polar 
coordinates (r, 9) in the dimensionless form are: 

1 T. -1 /^A U 2 \ 

Ut + UUr + - VUe = -Vr + Re Am o ve ] , 

ry iy \ ly^ ry ^ j 

1 uv 1 -^ ( K V 2 \ , , 

vt^uvr^ - vveA = — + Re Af -\ — uq \ , (1.1) 

ry ry ry \ cyii ry ^ j 

{ru)r + ve = 0, 

A( ) = ( )„. + - ( )r + \{ )ee. 

ry ry ^ 

Here u is the radial velocity, v is the azimuthal velocity, p is the pressure. Re is the Reynolds 
number. 

Boundary conditions at r = 1 and r = h have the form: 

w|.=i = 0, t^|,=i = l, = 0, vl^, = Q. (1.2) 

Dimensionless variables t, r, u, v, p correspond to dimensional variables t', r', u', v', p'\ 

t' = — , r' = rri, {u' ,v') = ilri{u,v), p' = ppii^rf, 

Re= — e = Re'\ b=—. (1.3) 
u ri 
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where u is the kinematic viscosity, ri is the radius of the inner circle, r2 is the radius of the outer 
circle, is the azimuthal velocity of the inner circle rotation, p is the constant fluid density. 
The problem p.ll) - (11.2p has the stationary axially symmetric solution (Couette flow): 

a 6^ 1 

u = 0, v = Voir), Voir) = -+Pr, « = P = (1-4) 

Let \l/(t,r, ^) be a stream function and ip{t,r,6) be a perturbation of the stream function: 

1 1 

U = --11)0 = --^!e, V = 'ipr + Vo = -^r- (1.5) 
If we substitute to ([ri|)-([Il2]), we obtain: 

A^i + ^(^,(A^),-^,(A^/.),) =£A2^-^(A^),, (1.6) 
V'L=i = 0, ^.|,=i = 0, V^|,=, = 0, V^.|,=, = 0. (1.7) 



2 Galerkin approximations 

To obtain the approximate solution of the problem f ll.6l) - p.7l) . Galerkin method is used. We 
approximate the stream function perturbation ip as follows: 

L N 

^ = 2 ^ ^ {{xnin{t) cos{n9) - Hmnit) sin(n6')} fmn{r). (2.1) 

m=l n=l 

Here, the functions fmn{r) satisfy the boundary conditions: 

/^„(1) = 0, /;„(!) =0, fmn{b)=0, /;„(&)= 0, (2.2) 

The standard Galerkin procedure reduces the problem p.GD - p.Tp to the system of the ordinary 
differential equations for Galerkin coefficients Xmnit), ymnif)- We also assume that the functions 
fmnij) satisfy the following conditions: 

6 

j fmn{r)Anfjn(,r)rdr = -6jm (2.3) 
1 

1 

A„( ) = (),, + -( ),--(). (2.4) 

In the work ^9j the two-dimensional flow in the annular strip (driven by the azimuthal force 
changing the sign with radius) is studied. Analogously to the work [19] the functions fmnij) for 
fixed n are taken as the eigenfunctions of the following boundary value problem: 

A2/ + /^2a„/ = 0, (2.5) 
/(1) = 0, /'(1) = 0, /(6) = 0, /'(6) = (2.6) 

Calculation of the eigenvalues hmn for the problem f l2.5p -f l2T6l ) and of the corresponding 
eigenfunctions fmn{r) is a tiresome procedure. That is why we also use the polynomial basis to 
construct functions fmn{^)- 
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To determine Galerkin coefficients Xmn{t), ymn{t), the following system of the ordinary 
differential equations is obtained: 



, L L L L N N 



j = l jr = l j = l p=l fc=l 5=1 



L L L L N N 



j=l j=l j=l p=l k=l q=l 

Here 

C(kqn = 1, if n = \k — q\ or n = k + q; auqn = 0, otherwise; 
Pkqn = 1, if n = |A; - g|; (3kqn = -1, if n = k + q; f3kqn = 0, otherwise. 
The coefficients 6mjn, Imjn and Kjkpqmn are defined by the following relations: 



6 6 

4„,„^//™.A„,,„r,, ,„,„^„/v„/„„,,„*. ,2.8, 



1 

6 



Kjkpqmn Q' j {.f jkdpq 9 jkfpi) fmn dv , (^jn ^nfjn- 



3 Numerical experiments 

Solutions of the system of equations (12.71) depend on the two dimensionless parameters: the 
Reynolds number Re (or e = Re~^) and on the ratio of the outer and inner radii b. Note that the 
parameter e is included to the system (12. 7p explicitly whereas the coefficients 6mjn, Imjn h Kjkpqmn 
depend on the parameter h. 

3.1 Basis functions 

Two sets of the basis functions are in use: the set of solutions of the problem ( 12.5p - (12.6p . and the 
linear combinations of the polynomials Pm{f^) = {r — 1)^(& — r^r"^'^ which satisfy the conditions 

o-dsai). 

Solutions of the problem (l2.5p -( l2T6ll have the form [T9| : 

fmnij'^ [^mn'-^n(p) ~l~ ^mn^ni^P) ~^ CmnP ~l~ DmnP '\ / ^mm P ^mn^ 

where J„ is the Bessel function of the first kind, F„ is the Neumann function. 

The eigenvalues h^n and coefficients A^n, -Bmn, Cmn, D^n are defined by the boundary 
conditions: 

fmn{l)=0, /;„(!) =0, fmn{b)=0, /^J^) = (3.1) 

Values hmn, (m = 1, 2, . . . , L) are the roots of the determinant of the system (13. ip . The nontrivial 
solution {Ajnn, Bmn, Cmni, Dmn) of the systcm (I3.ip must satisfy the condition ( 12.31) . The function 
Qmn has the form: 

9mn ^nfmn AmnJni^P^ ~^ Bmn^ni^P^ ■ 
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The system of functions Pm{f^) = {r — 1)^(6 — r)^r™~^ is used to construct the polynomial basis. 
Functions fmn{f) for fixed n are constructed with the help of Gramm-Shchmidt orthogonalization 
with the scalar product: 

{U,V) = -J^ u{r)AnV{r)rdr, {fmnjran) = 1- 

Chosen sets of basis functions are almost identical. Indeed, as the numerical experiment showed, 
properties of the solution of f ll.6l) - p.7p obtained by Galerkin method weakly depend on the choice 
of system of the basis functions fmnij'). The number L of the radial basis functions is chosen to be 
2 or 3, and the number of the azimuthal basis functions is not more then 9. 

3.2 Dependence of the critical Reynolds numbers on the number of basis 
functions 

Equations ( 12. 7^ have the solution which corresponds to the base regime (Couette flow). For the 
small Reynolds numbers this flow is stable, however, it can become unstable when the Reynolds 
number increases. The neutral curve Re = Reo(&) for = 9, L = 3 is shown in Fig. [1] (left). The 
neutral curve corresponds to the oscillatory instability with the azimuthal wave number n = 1. 
Neutral curves with the azimuthal wave numbers n = 2, 3 have similar form. 
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Figure 1: Left: the neutral curve Re = Reo(fe) of the Couette flow instability, corresponding to the 
azimuthal wave number n = 1 (L = 3, N = 9, polynomial basis). Right: dependence of the critical 
Reynolds numbers on the number of basis functions L at 6 = 5. Continuous line corresponds to 
computation with polynomial basis; circles correspond to basis composed of the eigenfunctions of 
(I2.5p - (l2.6p : azimuthal wave number n = 1. 

When the parameter b is fixed, the critical Reynolds numbers Reo(&) are strongly dependent on L 
(the number of radial basis functions). Numerical investigation of the linear stability of the solution 
Xmn = 0, Umn = of cquatious (12.71) shows that the critical Reynolds numbers are monotonously 
increasing with the increase of L when L ^ 6. Results of calculation at 6 = 5 are presented in Fig.[T] 
(right). Continuous line corresponds to the polynomial basis, and circles represent basis composed 
from the eigenfunctions of the problem (I2.5p - (l2.6p . 

3.3 Flow regimes 

Two different stable stationary regimes (traveling waves) are found apart from the base Couette 
flow in numerical experiments performed. Galerkin coefficients Xmn{t) and ymn{t) are periodic 
functions of time t, while the amplitudes dmn = {^mn + VmnY^'^ constant for these regimes. 
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One of them appears as a result of instability of the base Couette flow. Nature of the second regime 
is still not clear. Both flows have a number of typical features that are independent of the outer 
radius 6, the choice of the basis functions fmn{r) and the number of these functions L (for small 
enough L) provided that a number of the azimuthal basis functions is chosen to be large enough. 



0.4- 



0.2- 








Figure 2: Dependence of the amplitude dn on the Reynolds number: continuous line corresponds 
to regime with precessing vortex; circles correspond to regime with the jet. Left picture correspond 
to 6 = 5, right picture corresponds to 6 = 7, polynomial basis is used in both cases. Regime J is 
stable at 110 < Re < 290 for b = 7, and at 175 < Re < 480 for 6 = 5. 



In the wide interval of parameters regimes observed are stable. The onset of one or another of 
these regimes depends on the initial conditions. One regime is the stationary flow with the single 
precessing vortex spot. We call it 'regime V. Main feature of the second flow regime is the region of 
intensive flow in the form of the radial jet, precessing in the azimuthal direction. We call it 'regime 
J'. 

In the models constructed with the use of polynomial basis by r, the system loses stability 
smoothly in the sense that the periodical motion branches to the supercritical area Re > Reo and is 
stable. If we use the basis composed from the eigenfunctions of the problem f l2.5p - (12.6p . we observe 
the rough loss of stability: the unstable limiting cycle branches to the subcritical area. 
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Figure 3: Dependence of precession speed on the Reynolds number: upper points correspond 
to regime with precessing vortex, lower points correspond to regime with the jet. Left picture 
corresponds to 6 = 5, right picture corresponds to 6 = 7. Polynomial basis is used in both cases. 



We emphasize again that there is a certain range of parameters where both stationary flow 
regimes are stable. The onset of one or another flow regime depends on initial conditions. Regime 
V typically appears when the initial values of | v| are small. Regime J is usually observed when 
the initial values of |v| are large. Note that the values of dmn that correspond to regime J are 
significantly larger than ones for the fiow of V type. The amplitude curves dn — dniTi^G^ are 
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presented in Fig. El Note that regime V branches from the base Couette flow whereas regime J 
appears 'from the sky-blue'. When the Reynolds number increases, the oscillatory instability of 
flow regime J is observed. 

Another characteristic of the flow regimes V and J is the precession speed in the azimuthal 
direction w = uj/{2'k). Dependence of precession speed on the Reynolds number for two values of 
h is presented in Fig. [3l 



3.4 Flow with a precessing vortex (regime V) 

One of the flows observed in numerical experiments is the regime V: periodic flow of the type 
of a traveling wave with the single vortex that slowly precesses in the azimuthal direction (the 
direction of the vortex precession coincides with the direction of the inner circle rotation). Onset 
of this regime is linked to the loss of stability of the base Couette flow. Regime V appears for 
both sets of basis functions: a set of the eigenfunctions of the problem (12. 50 - 112. 6p and a set of the 
polynomials. Typical results of numerical experiment for the flow of the type V are presented in 
Figs. SHHl Computation is performed at L = 3, = 9, basis composed of the eigenfunctions of 
is used. 




Figure 4: Regime V, streamlines (1 
(3, 4); 1 and 3 correspond to 6 = 5, 



, 2) and level lines of the stream function perturbation ?/^(r, 6) 
Re = 320; 2 and 4 correspond to 6 = 7, Re = 200. 




Figure 5: Regime V, 6 = 5, Re = 320. Left: graph of the absolute velocity |v(r, 6*)!. Right: contour 
plot of the absolute velocity. 

The streamlines and the level lines of the stream function perturbation are plotted in Fig. [4] 
(location of the soliton-like vortex is clearly visible). Graph of the absolute velocity |v(r, 6*)! for 
the flow with the single vortex (regime V) is exposed in Fig. [5] (left). Contour plot of the absolute 
velocity in the rectangle r G 9 G [0, 27r] is shown in Fig. [5] (right). Regions of the different 

flow intensity are painted with different colors. The value |v| is taken as the measure of the flow 
intensity. The highest flow intensity is observed in the area close to the inner boundary (black 
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color). Note that the values of fluid velocity at the inner and outer boundaries are |v| = 1 and 
|v| = correspondingly, and |v| = k/7 {k = 1,2, ...,6) for level lines in Fig. 5 (right). The contour 
plots of radial velocity u{r,9) and azimuthal velocity v{r,9) are presented in Fig. [6l 




Figure 6: Contour plots of the radial and azimuthal velocity u{r,6)^ v{r,9). Flow with precessing 
vortex; 6 = 5, Re = 320 




Figure 7: Dependence of radial velocity u{r,9) and of azimuthal velocity v{r,9) on 6' at r = 1 + 
(3A; — 1)(6 — l)/20, where k is the number of the line. Flow with precessing vortex. 




Figure 8: Dependence of vorticity Aip{r, 9) and of stream function perturbation ip{r^ 9) on 9 at 
r = 1 + (3A; — l)(h — l)/20, where k is the number of the line. Flow with precessing vortex. 

In Figs. [THH] radial and azimuthal velocity profiles, vorticity profiles, and stream function 
perturbation profiles in the cross-sections r = const are presented. Note that zeros of the radial 
velocity u{r, 9) almost coincide with the extreme points of the azimuthal velocity f (r, 9), and zeros 
of A?/'(r, 9) almost coincide with zeros of ?/'(r, 9). Extreme points of A?/'(r, 9) and extreme points of 
ip{r, 9) also almost coincide. 



Low-order models of 2D Quid Row in annulus 



9 



3.5 Flow with a precessing jet (regime J) 

Apart from the regime with precessing vortex, another type of stationary flow is observed in 
numerical experiments. It is the slowly precessing in azimuthal direction stationary structure, whose 
main feature is existence of the jet-like area of intensive fluid motion. 

Results of numerical experiments for regime J are presented in Figs. [9HT41 The streamlines and 
the level lines of the stream function perturbation are plotted in Fig. [9l There is also the precessing 
vortex, but its shape is more complicated than that for the regime V. The vortex is stretched in 
the azimuthal direction. 




Figure 9: Regime J, streamlines (1, 2) and level lines of the stream function perturbation iplr, 9) (3, 
4); 1 and 3 correspond to & = 5, Re = 320; 2 and 4 correspond to 6 = 7, Re = 200. Basis composed 
of the eigenfunctions of f l2.5p -( l2T6l l is used in both cases. 



Graph of the absolute velocity |v(r, 6')| for the flow with precessing jet (regime J) and contour 
plot of the absolute velocity in the rectangle r e [1, 6], 6^ e [0, 27r] are shown in Fig. [10] (as before, 
|v| is taken as the measure of the flow intensity). The most intensive flow is observed in the certain 
region close to the inner boundary. Comparison with Fig. [5] shows that the boundary layer near 
the inner boundary have qualitatively different form for regimes V and J. 




Figure 10: Regime J, 6 = 5, Re = 320. Left: graph of the absolute velocity |v(r, 6')|. Right: contour 
plot of the absolute velocity. 



In Figs. fT0] - [T4l various features of the flow regime J are illustrated (h = 5, Re = 320). The 
contour plots of radial velocity u(r, 6) and azimuthal velocity f (r, 6) are presented in Fig. [TTl The 
jet-like area of intensive fluid motion is clearly visible in Fig. [TT] (left). 

Radial profiles of the absolute velocity |v(r, 9) \ for fixed r are plotted in Fig. [I2] (left) , azimuthal 
profiles of the absolute velocity for fixed 9 are presented in Fig. [I2] (right). Azimuthal profiles 
of velocities u{r,9) and v{r,9), profiles of vorticity Ailj{r,9) and profiles of the stream function 
perturbation ipi^r, 9) for r = const are plotted in Figs. [T3HT41 
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Figure 11: Contour plot of the radial velocity u{r,9) and of the azimuthal velocity v{r,9). Flow 
with precessing jet; 6 = 5, Re = 320. 

Profiles of azimuthal velocity for constant r chosen near to the inner or outer circle have the 
shock-wave type. The jump from larger to smaller values of azimuthal velocity in the direction 
of precession is observed near the inner circle while the jump from smaller to larger azimuthal 
velocities is observed near the outer circle. Profiles of radial velocity have the soliton-like shape 
for all fixed values of r. Note that zeros of the radial velocity u{r, 6) almost coincide with extreme 
points of the azimuthal velocity v(r, 6). 




Figure 12: Left: radial profiles of the absolute velocity |v(r, ^)| at ^ = (A; — l)vr/3, where k is the 
number of the line. Right: azimuthal profiles of |v(r, 9)\ at r = 1 + (S/c — l){h — l)/20, where k is 
the number of the line. Flow with precessing jet. 




Figure 13: Dependence of radial velocity u{r, 6) and of azimuthal velocity v{r, 9) on 9 at r — 
1 + {3k — l){b — l)/20, where k is the number of the line. Flow with precessing jet. 
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Figure 14: Dependence of vorticity Aip{r, 6) and of stream function perturbation ?/'(r, 6) on 6 at 
r = 1 + (3A; — l){h — l)/20, where k is the number of the line. Flow with precessing jet. 

It is interesting to compare results of computation with results of the experiments in which the 
jet-like flow regimes are observed. In Fig. [I5] (published with permission of V. A. Vladimirov and 
P. V. Denissenko), visualization of the experiment with the jet-like structure is presented (Re ~ 56). 
The qualitative similarity of the experimental results (Fig. [151 Isft) and of the numerical results 
(Fig. [TH right) is obvious. We repeat once again that the quantitative comparison can not be 
performed because the physical experiment is essentially three-dimensional while the computation 
considered in this paper is two-dimensional. 




Figure 15: Left: experiment. Free upper surface. Water-glycerol mixture, axial gap thickness 16 ■ 
lO^^m, radius of the cylinder 10~^ m, outer radius of domain 23.2 ■ 10"^ m, kinematic viscosity 
0.43-10~'^m^/s, rotation speed of the inner cylinder 230 rpm, Re ~ 56. Right: Regime J, streamlines, 
5 = 5, Re = 320. 
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